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The CKMT model for the nucleoli structure function F2 is in good agreement 
with the HERA data at low and moderate Q 2 . The fit to the same data obtained 
with a modified version of the model in which a logarithmic dependence on Q 2 
has been included is also discussed. Finally, we show how the parametrization of 
the CKMT model for the nucleon structure function F2 describes the HERA data 
when presented in the Caldwell-plot as the behavior of the logarithmic slopes of 
F2 vs x and Q 2 . 



1 The CKMT model 

The CKMT model for the parametrization of the nucleon structure function 
F2 is a theoretical model based on Regge theory which provides a consistent 
formulation of this function in the region of low Q 2 , and can be used as a safe 
and theoretically justified initial condition in the perturbative QCD evolution 
equation to obtain the structure function at larger values of Q 2 . 

The CKMT model El proposes for the nucleon structure functions 

F 2 (x, Q 2 ) = F s (x, Q 2 ) + F NS (x, Q 2 ), (1) 

the following parametrization of its two terms in the region of small and mod- 
erate Q 2 . For the singlet term, corresponding to the Pomeron contribution: 



F s (x, Q 2 ) = A- x'^ ■ (1 - *)"W 2 )+ 4 • 



1+A(Q 2 ) 

(2) 
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where the x^O behavior is determined by an effective intercept of the Pomeron, A, 
which takes into account Pomeron cuts and, therefore (and this is one of 
the main points of the model), it depends on Q 2 . This dependence was 
parametrized El as : 

Thus, for low values of Q 2 (large cuts), A is close to the effective value found 
from analysis of hadronic total cross-sections (A~0.08), while for high val- 
ues of Q 2 (small cuts), A takes the bare Pomeron value, A~0.2-0.25. The 
parametrization for the non-singlet term, which corresponds to the secondary 
reggeon (f, A2) contribution, is: 



F N s(x ) Q 2 ) = B-x l -^-(l-xr^^\^- l j , (4) 

where the x—>0 behavior is determined by the secondary reggeon intercept 
an, which is in the range afl=0.4-0.5. The valence quark contribution can 
be separated into the contribution of the u, B Ul and d, B e i, valence quarks, 
the normalization condition for valence quarks fixing both contributions at 
one given value of Q 2 (we use Qo = %-GeV 2 in our calculations). For both 
the singlet and the non-singlet terms, the behavior when x— >1 is controlled by 
n(Q 2 ), with n(Q 2 ) being 

^ = ? • fi + (5) 



2 

so that, for Q 2 =0, the valence quark distributions have the same power, given 
by Regge intercepts, as in Dual Parton ModelH, n(0)=an(0)— aAr(0)~ 3/2, 
and the behavior of n(Q 2 ) for large Q 2 is taken to coincide with dimensional 
counting rules. 

The total cross-section for real (Q 2 —0) photons can be obtained from the 
structure function F2 using the following relation: 

4TT 2 a E M „ ( „ 2 \ 



F 2 (x,Q 2 



(6) 

Q 2 =0 



The proper F2(x,Q 2 )~Q 2 behavior when Q 2 — >0, is fulfilled in the model due 
to the last factors in equations || and |J. Thus, the a^iy) has the following 
form in the CKMT model: 

a l ^{v) = A7T 2 a E M ■ (A ■ ■ (2m^) A ° + (B u + B d ) ■ b^* ■ {2mu) aR ' 1 ) . 

(7) 
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The parameters were determined^ from a joint fit of the <r*°* data and the 
NMC datacl on the proton structure function in the region lGeV 2 < Q 2 < 
5GeV 2 , and a very good description of the experimantal data available was 
obtained. 

2 Structure functions at high Q 2 

The next step in this approach is to introduce the QCD evolution in the par- 
tonic distributions of the CKMT model and thus to determine the structure 
functions at higher values of Q 2 . For this, the evolution equation in two loops 
in the MS scheme with A = 200MeV was usedtJ. 

As starting point for the QCD evolution, the value Q\ — 2GeV 2 , where 
the logarithmic derivative in Q 2 of -£2(2;, Q 2 ) in the model is very close to that 
obtained from the QCD evolution equation, was chosen. 

The|-results obtained by taking into account the QCD evolution in this 
way are 111 in a very good agreement with the experimental data on F2(x, Q 2 ) 
at high values of Q 2 . 

Although the fit of the NMC data was restricted to the region lGeV 2 < 
Q 2 < 5GeV 2 , it is interesting to mention that a good fit can also be obtained 
with the model up to Q 2 ~ lOGeV 2 . This allows to start the perturbative 
QCD evolution at larger values of Q 2 . 

3 Description of the HERA data on F 2 at low Q 2 

When the CKMT model was first usedEI to fit the available experimental data 
on i*2, the lack of experimental data at low and moderate Q 2 limited the 
accuracy in the determination of the values of the parameters in the model. 
Later on, the publication of the new experimental datatm on F2 from HERA 
at low and moderate Q 2 provided the opportunity to include in the fit of 
the parameters of the model experimental points from the kinematical region 
where the paramctrization should give a good description without need of any 
perturbative QCD evolution. n 

Thus, one proceededu as one had done in the previous fitEl, but by adding 
the above mentioned experimental data on F2 from HI and ZEUS at low^md 
moderate Q 2 , to thosetl from NMC and E665 collaborations, and to dataLj on 
cross-sections for real photoproduction, into a global fit which allowed the test 
of the model in wider regions of x and Q 2 . One took as initial condition for 
the values of the different parameters those obtained in the previous fit El. The 
result of the new common fit to <r'°* and F2 is presented in figures [l] and |[ 
and the final values of the parameters can be found in Table^(b) . The quality 
of the description of all the experimental data provided by the CKMT model, 
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and, in particular, of the the new experimental data from HERA is very good, 
with a value of x 2 /d.o.f. for the global fit, x 2 /d.o.f. =106. 95/167, where the 
statistical and systematic errors have been treated in quadrature, and where 
the relative normalization among all the experimental data sets has been taken 
equal to 1. Moreover, since the small- a; HERA experiments allowed for the first 
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Figure 1: u^ al and it*°*° 1 in ^tbarns vs W in GeV for different values of Q 2 . Theoretical 
fits have been obtained with the CKMT model (full line) and the modified version of the 
CKMT m odel (dashed line). Points at (a), Q 2 =0.GeV 2 (*8.), (b), Q 2 =0.15GeV 2 (*6.), (c), 
Q 2 =0.25GeV 2 (*5.), (d), Q 2 =0.5GeV 2 (*4. ), (e), Q 2 =0.8GeV 2 (*3.), (f), Q 2 =l.5GeV 2 
(*2.) , and (g), Q 2 =3.5GeV 2 (*1.)- Experimental data on F% arc from references 

(3), (4), and experimental data on a^ al are from references (7). 

time the experimental study of the question of the interplay between soft and 
hard physics, the model, which basically has only power dependence on Q 2 , 
was modified to include a logarithmic dependence on Q 2 as the one predicted 
asymptotically by perturbative QCD. 

To include the logarithmic dependence on Q 2 in our model, we take into 
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account that the behavior of F2 at small- a: in QCD is given by the singularities 
of the moments of the structure functions u, the .rightmost singularity giving 
the leading behavior. Thus, the following factors which correspond to the 
moments of the structure functions in the language of the OPE expansion, 
and can be calculated by the convolution in rapidity of the hard-upper part 
with the soft-lower part of the leptoproduction diagram, were introduced in 
the expression that the CKMT model gives for F^: 



«.(Qo) 



i = S, NS, (8) 



a s {Q 2 ) 

where the strong coupling constant is taken as 

Air 

<*s(Q) = 7 r, (9) 

\ Ji QCn J 

with M~lGeV, a hadronic mass included in ^| to avoid the singularity in a s 
when Q 2 — > Aq C£) , Aqcd=0.2 GeV, and (3o=ll— |n/ (in the calculations a 
number of flavors "./=3 was used). The exponents ds(n>s) and d,Ns( n Ns) m H 
are proportional to the largest eigenvalue of the anomalous dimension matrix, 
and to the anomalous dimension, respectively: 

Mns)~^— } - d , (10) 



with 

"""ft'" 1- A) 



do = f,d 1 = 1 -^l, (11) 



and 

4sKs) = ^ — t, (o 7 r- T\ + 7 ~ s i( n Ns) ) , (12) 

33 -2n f \2n N s{n N s + 1) 4 



with 



Si {n NS ) = n NS ■ V — — -. (13) 

^ k{k + n NS ) 



k=l 



Thus, we modify [l] in the following wayi 



(14) 
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Figure 2: F2(x,Q 2 ) vs Q 2 (in GeV 2 ) for different values of x. Theoretical fits have been 
obtained with the CKMT model (full line) and the modified version of the CKMT model 
(dashed line). Experimental points at (a), from left to right, x=0.4210~ 5 , x=0.44-10 — 5 , 
and x=0.46-lCr 5 (*8.); (b), from left to right, x=0.85-10~ 5 , x=0.8410~ 5 , x=0.831CT 5 , 
and x=0.86-10~ 5 (*6.) ; (c), from left to right, x=0.13-lCT 4 , and thr ee points at 
x=0.14-lCT 4 (*5.); (d), x=0.5-10 -4 (*4.); (e), x=0.8-lCr 4 (*3.); (f), x=0.2-10" 3 (*2.); (g), 
x=0.5-10~ 3 (*1.); (h) x=0.52-10" 2 (*1.); (i) x=0.13-10- 1 Experimental points for F 2 

are from references (3), (black circles), (4), (crosses), and (6), (black squares for the E665 
points, and black diamonds for the NMC points). 

The exponents of the new factors in [li], ds(ns) and djvs(^ivs); gi ye us the 
singularities in rii, i=S,NS, of the momenta, which, as we mentioned above, 
control the QCD small- a: behavior of F2. Therefore, in the CKMT model, 
these exponents have to be evaluated (see || and ||), at ?i5=l+A(Q 2 ^oo)= 
l+Ao(l+Ai), and at n^s—^R.-, respectively. 

Then, this modified version of the CKMT parametrization of F2 was used 
to repeat the fit of the same experimental data, including the HERA data 
on F2 at small and moderate Q 2 . As starting point for the QCD evolution, 
one takes the same value that Q\=2.GeV 2 was used to fix the normalization 
of the valence component. The result of this second fit is also presented in 
figures and 0, and the final values of the parameters in the model are given 
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in TableQ(c). 

As it can be seen in the figures [l] and |[ the quality of this second fit 
is also reasonable, although the value of x 2 /d.o.f. (x 2 /d.o./.=453. 19/167), is 
now appreciably higher than in the fit obtained with the non-modified version 
of the CKMT model. 



Table 1: Values of the parameters in the CKMT model obtained in former fits, (a), in the fit 
in which also the low Q 2 HERA data have been included, (b), and in the fit to the same data 
obtained with the modified version of the CKMT model in which a logarithmic dependence 
of F2 on Q 2 has been taken into account, (c). All dimensional parameters are given in GeV 2 . 
The valence counting rules provide the following values of B u and Bj, for the proton case, 
when fixing their normalization at Q 2 =2.GeV 2 : (a) B„=1.2064, B d =0.1798; (b) B n =1.1555, 
B d =0.1722; (c) B u =0.6862, B d =0.09742. In previous fits, (a), the parameter Ai had been 
fixed to a value Ai=2. 



CKMT model 


(a) 


(b) 


(c) 


A 


Q.1502 


0.1301 


0.1188 


a 


0.263f 


0.2628 


0.07939 


A 


0.07684 


0.09663 


0.1019 


Ai 


2.0 (fixed) 


1.9533 


1.2527 


A 2 


1.1170 


1.1606 


0.1258 


c 


3.5489 


3.5489 (fixed) 


3.5489 (fixed) 


b 


0.6452 


0.3840 


0.3194 


a R 


0.4150 


0.4150 (fixed) 


0.5872 



4 Description of the Caldwell-plot 

The so-called Caldwell-plot shows the logarithmic slope of the structure 
function F2, dF2/dlnQ 2 , derived from the ZEUS data, as a function of x, by 
fitting F-2 ~ a + blnQ 2 in bins of fixed x, using only statistical errors. This 
plot allows the study of the QCD scaling violations of F2, and, in particular, 
in the small- a: domain now accessible at HERA, where dF2/dlnQ 2 is directly 
related to the gluon density, can be an useful tool to investigate down to 
which value of Q 2 the perturbative NLO DGLAP QCD predictions give a 
good description of the F2 data. Thus, the logarithmic slope dF2/dlnQ 2 can 
be used to investigate the fundamental question of the interplay between the 
soft and the hard physics. 
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As it can be seen in reference (10), for values of x down to 3 • 10~ 4 , the 
slopes are increasing as x decreases, but at lower values of x and Q 2 the slopes 
decrease, what seems to indicate a deviation from the perturbative behavior 
of the hard regime. 

Also in the reference0 one can see how both Regge based parametrizations 
with a constant effective value for the Pomeron intercept, like the Donnachie- 
Landshoff Regge fit or the ZEUSREGGE fit, and pure perturbative NLO QCD 
predictions, like GRV94 NLO QCD fit or the ZEUSQCD fit, fail in describing 
correctly the experimental data in the whole kinematical region. While the 
Donnachie-Landshoff and the ZEUSREGGE fits do not describe the data for 
values of x larger than ~ 10~ 5 , GRV and ZEUSQCD do not follow the data 
when one goes to values of x smaller than ~ 6 ■ 10~ 5 (see (10) and references 
therein for more details on this discussion). 

On the other side, the CKMT model described above, based on_the Regge 
behavior, but with a Q 2 -dependent Pomeron intercept, describes 113 the data 
in the region of low Q 2 , and when taken as the initial condition for the QCD 
evolution equations, provides a complete description of the experimental results 
in the whole ranges of x and Q 2 . 

As a matter of fact, by using the formulae of the pure CKMT model in 
section 1, when x is kept fixed one can writeO for the CKMT model the slope 
dF 2 /dlnQ 2 as: 



'-^if 1 = Fs(x, Q 2 )[q^k- 2 (A(Q 2 ) - Ao) 



- (n(Q 2 ) - |) In(l -x) + 



F NS (x,Q 2 )[Q^- c (n(Q 2 )-l)ln(l-x) (15) 

Q 2 +b J' 



that in the limit Q 2 — > takes the form 



d J^~ { l + A )F s ( x ,Q 2 ) 

+ a R (0)F NS (x, Q 2 ). (16) 



Also, if one considers the case when W is fixed one can take x ~ cte ■ Q 2 , and 
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then, up to constant factors, one gets: 



^%fl = F s (x, Q 2 )[-^fe (A(Q 2 ) - A ) HQ 2 + a) 
-A(Q 2 ) + -^r c (n(Q 2 ) - §) ln(l - Q 2 ) 

Q 2 n(Q 2 ) a(l+A(Q 2 )) 1 
1-Q 2 ~ l ~ Q 2 +a J (17) 

+ F^(ar, Q2)[_^_(„(Q2) - I )l„(i - Q*) 

+S + (i-M0))- 2 S 1 ]. 



Now, if one takes W fixed with Q 2 ~ x — > 0, one can easily see that this 
equation simply reduces to: 



dF 2 {x,Q 2 ) 
dlnQ 2 



F 2 (x,Q 2 ). (18) 



One has to note that both equations 16 and |l8| are valid for any well-behaved 
parametrization of F 2 (i.e., any parametrization fulfilling the relation in equa- 
tion ||) . 

Taking intxL account the general features of the CKMT model described 
above, we uselij the CKMT model to describe the experimental data in the 
region of low Q 2 (0 < Q 2 < Q 2 , = 2.GeV 2 ), and then we take this parametriza- 
tion as the initial condition at Qq = 2.GeV 2 , to be used in the QCD evolution 
equation to obtain a description of the experiment at values of Q 2 higher than 
Qq = 2.GeV 2 . We present our results in the shape of both the dF 2 / dlnQ 2 and 
the dlnF2 / dln{l / x) slopes in order to compare with the experimental data 
when presented in the so-called Caldwell-plot. 

The way we proceed to calculate F 2 , and the derivatives dF 2 /dlnQ 2 and 
dlnF2 j dln(l / x) is the following (see reference (12) and the appendix there for 
all the technical details on how the QCD evolution has been performed): 

• In the region < Q 2 < Q\ — 2.GeV 2 we use the pure CKMT model for 
F 2 . 

• For Qq < Q 2 < charm threshold^, we make the QCD evolution of F 2 
at NLO in the MS scheme for a number of flavours rif = 3, and we take 
as the starting parametrization the CKMT one at Qq = 2.GeV 2 . 



When charm threshold < Q 2 < Q 2 = 50.GeV 2 , also the QCD evolution 
of F 2 is implemented at NLO in the MS scheme for a number of flavours 
n/ = 3, using the parton distribution functions for the it, d, s quarks, and 
by including the charm contribution via photon-gluon fusion. 
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• For values of Q 2 > Q 2 , QCD evolution is computed at NLO in the MS 
scheme, but now with a number of flavours nj = 4, and by using the 
parton distribution functions for the it, d, s, and c quarks. 

One has to note that in the treatment of the charm contribution we have 
followed reference (13). 

The results we have obtained are presented in figures || to ||. In figure || 
(Caldwell-plot), the slope dF 2 /dlnQ 2 is shown as a function of x, and compared 
with the a + blnQ 2 fit to the ZEUS F 2 data in bins of x. 

dF 2 /dln(Q T ) 

0.5 I 1 r-n— , 1 — q 1 1 — q 1 ,-n— , r-T— . r-r. 



0.45 - 




-0.05 — J— ■ — J—- — J— ■ — J— ■ —J 

le-06 le - 0.5 0.0001 0.001 0.01 0.1 1 



Figure 3: dF2/dlnQ 2 as a function of x computed in the CKMT model (see reference (12) 
for details on the calculation), compared with the fit of the ZEUS Fi data in bins of x to 
the form a + blnQ 2 (see reference (10) and references therein for more details on the data 
and the experimental fit). 

Figures ||and|| show the slope dlnF2 / dln{l / x) as a function of Q 2 com- 
pared to the fits F 2 = Ae~ a<,// of the the ZEUS and HI data, respectively. In 
Figure |J, as the x range of the BPC95 data is restricted, also the E665Q data 
were included in (10), and are now also taken into account. The interest of 
these figures is clear, since this slope can be interpreted as the effective A of 
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the Pomeron exchange, A e // = dlnFI j dln(\ / x) . In the experimental fits, each 
Q 2 bin corresponds to a average value of a;, < x >, calculated from the mean 
value of ln(l/x) weighted by the statistical errors of the corresponding Fi val- 
ues in that bin. Even though we can proceed as in the experimental fits, and 
we get a very good agreement with the data, since the estimation of < x > is 
in some sense artificial and arbitrary, and introduces unphysical wiggles when 
drawing one full line connecting the different bins, we prefered to make for all 
the Q 2 bins in this figures the choice of the smallest x in the data instead of 
considering a different < x > for each Q 2 . This choice is based on the fact that 
the ansatz A e // = dlnF2 / dln(l / x) is actually valid for small x, and results in 
a smooth curve except for the jump in the region around Q 2 ~ 50GeV 2 , where 
the evolution procedure changes (again, see reference (12) for more details). 

dlnF 2 /dln(l/x) 

0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 


-0.1 

0.1 1 10 100 1000 

Q 2 

Figure 4: dlnFi j 'dln(l / 'x) as a function of Q 2 calculated in the CKMT model, and compared 
to the fit F 2 = Ax~ x "ff of the ZEUS (10) and E665 (6) data with x < 0.01 For details on 
the CKMT calculation, see reference (12). 

Finally, figure |^ is the compilation of the behavior of F2 as a function of 
x for twelve different values of Q 2 (from Q 2 = 0.6GeV 2 to Q 2 = U.GeV 2 ), 
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Figure 5: dlnF2 / dln(l / x) as a function of Q 2 calculated in the CKMT model, and compared 
to the fit F 2 = Ax~ X 'SS of the HI data (11). For details on the CKMT calculation, sec 
reference (12). 



corresponding to the values presented by the ZEUS Collaboration in reference 
(10). 

A very good agreement with the experiment is obtained for all the x and 
Q 2 values, showing that the experimental data can be described by using as 
initial condition for the QCD evolution equation a model of F% where the 
shadowing effects which are important at low values of Q 2 are included. 

5 Conclusions 

The CKMT model for the parametrization of the nucleon structure functions 
provides a very good description of all the available experimental data on 
F2(x,Q 2 ) at low and moderate Q 2 , including the more recent small- 2: HERA 
points. A second fit to the same data obtained with a modified version of the 
model in which a logarithmic dependence on Q 2 is included, has been also 
presented. Even though the quality of this second description is reasonable, 
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its x 2 /d.o.f. is appreciably higher than that corresponding to the fit obtained 
with the non-modified version of the CKMT model. 

Finally, the CKMT model for F-i has been used as the initial condition 
in the QCD-evolution equation to describe the HERA experimental data pre- 
sented in the so-called Cald well-plot, where the ^-dependence of the logarith- 
mic slope of the structure function, dF 2 /dlnQ 2 , is shown for different Q 2 bins, 
and in the plot of the Q 2 -dependence of the X e ff, i.e., of the Q 2 behavior of the 
slope dlnF%/dln{l/x), now for different bins of x. The obtained results show 
that the available experimental data can be described by performing the QCD 
evolution of a model of F2 where the shadowing effects which are important 
at low values of Q 2 are included. 
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Figure 6: F2 as a function of x computed in the CKMT model (12) for twelve different values 
of Q 2 , and compared with the following experimental data (see (10) for the experimental 
references): ZEUS SVX95 (black circles), HI SVX95 (white triangles), ZEUS BPC95 (white 
squares), E665 (white diamonds), and ZEUS 94 (white circles). The dotted line is the 
theoretical result obtained with the pure CKMT model, and the solid line is the result 
obtained with the QCD-evoluted CKMT model. 
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